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and neutrino energies for the oscillation probabilities to exhibit certain desirable features. 
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1 Introduction 



When performing long-baseline neutrino oscillation experiments on the Earth with ac- 
celerator based beams, or when detecting atmospheric neutrinos coming from below, the 
neutrinos necessarily traverse the Earth's interior [1-7]. This makes the understanding 
of matter effects [8-10] on the neutrino oscillation probabilities an indispensable part of 
analyzing such experiments. 

Though the exact three-flavor neutrino oscillation probabilities in constant-density 
matter can be expressed analytically, 1 the exact formula is too complicated to yield physical 
insight [12-15]. Consequently, various analytical approximations have been devised to 
better understand the physics potential of various neutrino experiments [9, 16-21]. These 
approximations relied on expansions in the small parameters a = Sm^i/dm^ ~ 0.03 
and/or sx3 = sin #13 in one form or another, a systematic treatment of which can be found 
in Ref. [19]. Of the formulae available, the one derived by Cervera et al. in Ref. [22] has 
been used frequently due to its simplicity. Unfortunately, the accuracy of the Cervera et 
al. formula suffers when the value of #13 is as large as was measured by Daya Bay [23] and 
RENO [24], consistent with both earlier and later determinations by T2K [25], MINOS 
[26, 27], and Double Chooz [28, 29]. This is due to the fact that the Cervera et al. formula 
includes terms of order 0(a 2 ), 0(as\3), and 0(sf 3 ). Given that the current world average 
of S13 = sin 613 is about 0.15 [30], not all terms of the same order are included. Asano and 
Minakata [31] have derived the order 0(asf 3 ) and 0(sf 3 ) corrections to the Cervera et al. 
formula, but the simplicity of the original expressions is lost. 

In previous papers [32, 33], we had advocated the utility of the Jacobi method [34] in 
deriving approximate formulae for neutrino oscillation in matter. The formula presented 
there maintained the expressions for the oscillation probabilities in vacuum, and only re- 
placed the mixings angles and the mass-squared differences with their effective values in 
matter. The approximation worked very well except when #13 was very small, a possibil- 
ity that could not be ignored until the Daya Bay/RENO measurements. In this paper, 
we reintroduce the method with further refinements which improve the accuracy of the 
approximation, while maintaining its ease of use. 

This paper is organized as follows. In section 2, we explain our approach to the matter 
effect problem, and list all the formulae necessary to calculate our approximation. In 
section 3, we demonstrate the accuracy of our approximation at various baseline lengths, 
different mass hierarchies, and different values of the CP violating phase 5. Comparisons 
with the approximations of Cervera et al. [22] and Asano-Minakata [31] are also made. In 
section 4, we show how simple calculations using our approximation can be used to derive 
the baselines and energies at which the oscillation probabilities exhibit desirable features. 
We conclude in Section 5. Detailed derivation of our approximation is given in appendices 
B and C. 



1 The cubic characteristic equation for the eigenvalues of the effective Hamiltonian can be solved analyt- 
ically using Cardano's formula [11]. 
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Figure 1. The dependence of the line-averaged mass density p on the baseline length L based on 
the Preliminary Reference Earth Model [35] . The labels on the right edge of the frame indicate the 
corresponding values of a/E. The green and red dashed lines indicate pL = 54000 km • g/cm 3 and 
pL = 32000 km • g/cm 3 , respectively, which are conditions that will be discussed in section 4. 

2 The Approximation 

In the following, we use the conventions and notation reviewed in Appendix A. 
2.1 Diagonalization of the Effective Hamiltonian 

If the matter density along the baseline is constant, the effective Hamiltonian which governs 
the evolution of neutrino flavor in matter is given by 

(2.1) 

where U is the neutrino mixing matrix in vacuum, and 

a = 2V~2G F N e E = 7.63 x 1CT 5 (eV 2 ) (-^ (JL^j . (2.2) 

Here, N e is the electron number density, p the matter mass density along the baseline, 
and E the neutrino energy. The above term appearing in the ee-component of H a is 
due to the interaction of v e with the electrons in matter via VF-exchange, and Eq. (2.2) 
assumes iV e = N p = N n in Earth matter. It also assumes E <C Myy since the VT-exchange 
interaction is approximated by a point-like four-fermion interaction. Z-exchange effects are 
flavor universal and only contribute a term proportional to the unit matrix to H a , which 
can be dropped. 
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If we write the eigenvalues of H a as A, 
U, that is 



1,2,3) and the diagonalization matrix as 



U 



Ai 
A 2 
A 3 



u 



(2.3) 



then the neutrino oscillation probabilities in matter are obtained by simply taking their 
expressions in vacuum and replacing the elements of the mixing matrix U and the mass- 
square differences Smfj with their "effective" values in matter: 



U, 



5ml 



6Xi 



Xi — A, 



(2.4) 



Note that a is -E-dependent, which means that both U a i and SXij are also E-dependent. 
They also depend on the baseline length L since the average matter density p along a 
baseline varies with L. The L-dependence of the average p and the corresponding value of 
a/E are shown in Fig. 1. 

For anti- neutrino beams, the flavor-evolution Hamiltonian in matter is 



(2.5) 



In comparison to Eq. (2.1), the CP violating phase 5 in U and the matter-effect term a 
both acquire minus signs. Let us write the eigenvalues of H a as Aj (i = 1,2,3) and the 
diagonalization matrix as U, that is 
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H a = U* 


5m\ x 


U T + 
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A 2 
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(2.6) 



Note that the tilde above U here is flipped to distinguish it from U in Eq. (2.3). The anti- 
neutrino oscillation probabilities in matter are then obtained by making the replacements 



5 m: 



5Xi 



Xi - X 



3 > 



(2.7) 



in the vacuum expressions. 

2.2 Effective Mixing Angles 

As mentioned in the introduction, it is possible to write down exact analytical expressions 
for U a i and SXij, as we h as their anti-neutrino counterparts [12]. However, simpler and more 
transparent approximate expressions can be obtained using the Jacobi method [34]. In our 
approach, the matrix U is approximated by taking the conventional parameterization of 
the matrix U in vacuum in terms of three mixing angles and a CP- violating phase, namely 



U = #23(023, 0)i?l 3 (013,<5)#12(012,O) , 



(2i 
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and replacing the values of 612 and #13 with their effective values in matter: 2 

U ~ R23(023,O)R 13 (e' 13l 8)R l2 (d' 12 ,O) . (2.9) 

The angle 623 and the CP-violating phase 5 are kept at their vacuum values. The eigen- 
values Xi (i = 1,2,3) of H a are also given approximate expressions. Thus, to take matter 
effects into account when calculating neutrino oscillation probabilities, all that is necessary 
is to take their expressions in terms of the mixing angles and CP-phase in vacuum as is, and 
replace the two angles as well as the mass-squared differences with their effective values in 
matter: 6*12 — > 0' 12 , #13 — > #13, — > 8\j = \ — ^j- This simplifies the calculation con- 
siderably, and allows for a transparent understanding of how matter-effects affect neutrino 
oscillation by looking at the o-dependence of the effective parameters. 
Similarly, in the anti-neutrino case, we approximate 

U « i? 23 (023, 0)R l3 (e' 13 ,5)R 12 (rf 12 , 0) . (2.10) 

Again, 823 and 5 are unaffected while O12 and #13 are replaced by their effective values in 
matter. The eigenvalues Aj (i = 1, 2, 3) of H a are given approximate expressions as in the 
neutrino case. Thus, the calculation of matter effects for anti-neutrino beams entails the 



-/ 



replacements 6*12 — > 12 , #13 — > 13 , Smfj — > 5Xij = Xi — Xj. 
2.3 Definition of (3 

In the following, we list the expressions necessary to calculate the effective mixing angles 
and the effective mass-squared differences for the neutrino and anti-neutrino cases sepa- 
rately. Detailed derivation of our approximation is given in Appendix B. We will also show 
plots depicting how various effective parameters depend on the matter-effect parameter 
a. Due to the wide separation in scale between bm^ x and Sm 31 , we find it convenient to 
introduce the parameter (5 via 3 



31 1 



l^ m 3ll V \^ m . 



and plot our effective parameters as functions of (3 instead of a. f3 = corresponds to 
a = l^mfj-J, j3 = — 2 to a = dm^i, and so on. 

2.4 Neutrino Case 

For the neutrino case, the effective mixing angles in matter are given by 

,2 /„2 



0Sm^/cj 3 ) sin 2fli2 
(om^/c^g) cos Z&12 — a 

tan 2^ = (^1 ; ) sin2 ^ 2) 

(oto^ — om^sy cos 20i3 — a> 



Absorbing matter effects into 812 and #13 is similar in spirit to Ref. [18]. 

3 We avoid the use of the symbols a or A since they often respectively denote Sm^i/Sm^ and a/Sm^ 
in the literature. 
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-4-3-2-1 1 2" "-4 -3-2-10 1 2 

P P 
(a) neutrino mixing angles (b) anti-neutrino mixing angles 

Figure 2. The dependences of the effective mixing angles on f3 = — log 6 (a/|Jm| 1 1) for the neutrino 
(a) and antineutrino (b) cases. (3 — corresponds to a — \dm 31 \, and f3 — —2 to a — Srn^i- The f3- 
dependences of 0' 13 and 9 13 depend on the mass hierarchy: when Sm 31 > (normal hierarchy, NH) 
9' 13 increases toward 7r/2 whereas 8 13 decreases toward zero, while in the Sm 31 < case (inverted 
hierarchy, IH), it is the other way around. 




P P 
(a) neutrino case (b) anti-neutrino case 

Figure 3. The /3-dependences of the sines of twice the effective mixing angles for the neutrino 
(a) and antineutrino (b) cases. The difference in the behavior of the effective #13 mixing angle for 
normal and inverted hierarchies will allow us to determine which is chosen by nature. 



The /3-dependence of these angles are shown in Fig. 2(a), and that of the sines of twice 
these angles in Fig. 3(a). 

The eigenvalues of the effective Hamiltonian are approximated by 

\i « A'_ , 
A 2 w A" , 

A3 « X'± , (2.13) 
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(a) neutrino, normal hierarchy 




P 

(c) anti-neutrino, normal hierarchy 



-4 -3-2-10 1 2 

P 

(b) neutrino, inverted hierarchy 




P 



(d) anti-neutrino, inverted hierarchy 



Figure 4. Dependence of the approximate eigenvalues of the effective Hamiltonian on j3 = 
— log e (a. / 1 (5tti§ 1 I ) for the (a) neutrino normal hierarchy, (b) neutrino inverted hierarchy, (c) anti- 
neutrino normal hierarchy, and (d) anti-neutrino inverted hierarchy cases. 



where the upper(lower) sign is for the normal (inverted) hierarchy, with 



A' 



± 



{8m 2 2l + ac 2 l3 ) ± y/{5- 



in 



21 



ac 2 l3 ) 2 + 4ac? 3 sf 2 <5mi 1 



4 = 



[A' + + (Sm^ + «s? 3 )] ± V [A'+ - (Sm 2 31 + asj 3 )] + 4a 2 s ,2 2 c 2 13 s 



13 



(2.14) 



and s' 2 2 = sin 2 8' 12 ■ The /3-dependence of A'_ and A " for the normal and inverted mass 
hierarchies are shown in Figs. 4(a) and (b), respectively. That of A^_ is shown in Fig. 12(b). 
For the inverted hierarchy case, 5m 31 < 0, the above expressions simplify to 

A 2 X" « A' + , A 3 ~ X'L « 5m 2 31 < . (2.15) 
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2.5 Anti-Neutrino Case 



For the anti-neutrino case, we need to flip the sign in front of the matter effect factor a, 
and the effective mixing angles are given by 



tan 29 12 



(5m| 1 /cf 3 ) sin 2^12 
(<frrc|i/cf 3 ) cos 2#i2 + a 
(5m 31 — Sm^s^) sin26> 



tan2#' 13 = rpi^ -—-13 (2 lg) 

(tajj — Om^sjj) cos 26^13 + a 

The /3-dependence of these angles are shown in Fig. 2(b), and that of the sines of twice 
these angles in Fig. 3(b). 

The three eigenvalues of the effective Hamiltonian are approximated by 

Ai ~ A T , 
A 2 ~ A + , 

A 3 « Al , (2.17) 
where the upper(lower) sign is for the normal (inverted) hierarchy, with 
{5m 2 21 - ac 2 l3 ) ± v 7 {5m 2 21 + ac 2 3 ) 2 - Aac{ 3 s 2 l2 5m 2 21 



2 



_ [\'_ + (5m 2 31 - as 2 13 )} ± ^ [X_ - ((Sm^ - as 2 3 )] 2 + 4a 2 c ,2 2 c 2 3 s 2 3 
A± = ^ , (2.18) 

and c'^ = cos 2 # 12 . The /3-dependence of A , and A± for the normal and inverted mass 
hierarchies are shown in Figs. 4(c) and (d), respectively. That of A_ is shown in Fig. 15(b). 
For the normal hierarchy case, Sm 31 > 0, the above expressions simplify to 

Ai « A" « A_ , A 3 « A" « 5m 2 3l . (2.19) 

2.6 Comment on Coding 

As seen above, our formulae for both the neutrino and anti-neutrino cases are fairly compact 
and easy to code. In particular, the effective mixing angles for the neutrino and anti- 
neutrino cases can be calculated with the same code by simply flipping the sign of the 
matter-effect parameter a, cf. Eqs. (2.12) and (2.16). The same can be said of A± and 
A ± defined in Eqs. (2.14) and (2.18). In the case of A± and A±, one also needs to make 
the swap A^ -R- A_ but otherwise the code will be essentially the same. Once the effective 
mixing angles and the A's are obtained, one only needs to input them into the code for 
calculating the vacuum oscillation probabilities, and we are done. 



3 Demonstration of the Accuracy of the Approximation 

In this section, we demonstrate the accuracy of our approximation. For the vacuum values 
of the mixing angles and mass-squared differences, we use the global fit values from Ref. [30] 
listed in Table 1. 



- 8 - 




Figure 5. Comparison of the approximation formulae of Cervera et al., Asano-Minakata, and this 
work at L = 4000km. The average matter density for this baseline length is p — 3.81 g/cm 3 . The 
dashed line gives the exact numerical result. [35] 



Neutrino, L-810km, (5-0, Normal Hierarchy 




0.002 



Neutrino, L-8I0km, 6=0, Normal Hierarchy 



Cervera et al. 
Asano-Minakata 
This work 




(a) 



(b) 



Figure 6. Comparison of the approximation formulae of Cervera et al., Asano-Minakata, and this 
work at L — 810 km, which is the distance from Fermilab to NO^A. The average matter density 
for this baseline length is p — 2.80 g/cm 3 . The dashed line gives the exact numerical result. [35] 



We begin by comparing our approximation to those of Cervera et al. [22] and of 
Asano and Minakata [31]. In Fig. 5(a), we plot the approximate — > u e oscillation 
probabilities calculated using the three approximations against the exact numerical result 
for the baseline length L = 4000 km. This is the distance used by Asano and Minakata in 





7.5 x 10" 5 eV 2 




2.47 x 10~ 3 eV 2 


sin 2 6>23 


0.5 


sin 2 6>i2 


0.3 


sin 2 0i3 


0.023 



Table 1. Best-fit values of oscillation parameters taken from Ref. [30]. 
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Figure 7. Comparison of the approximation formulae of Cervera et al., Asano-Minakata, and this 
work at L = 8770 km, which is the distance from CERN to Kamioka. The average matter density 
for this baseline length is p = 4.33 g/cm 3 . The dashed line gives the exact numerical result. Note 
that the Asano-Minakata formula gives negative probability for E ~ 4 GeV. [35] 



Neutrino, L= 1300km, Normal Hierarchy 



Neutrino, L= 1300km, Inverted Hierarchy 




E (GeV) 

(a) 

Anti-neutrino, L=1300km, Normal Hierarchy 





Anti-neutrino, Z,= 1300km, Inverted Hierarchy 




Figure 8. Comparison of our approximation formulae (colored) to the exact numerical results 
(black, dashed) for various values of the CP violating phase S at L = 1300 km. The average matter 
density at this baseline length is p = 2.87 g/cm 3 . 
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Neutrino, L=2300km, Normal Hierarchy 



Neutrino, L=2300km, Inverted Hierarchy 





(a) 



(b) 



Anti -neutrino, L=2300km, Normal Hierarchy 



Anti -neutrino, L=2300km, Inverted Hierarchy 





(d) 



Figure 9. Comparison of our approximation formulae (colored) to the exact numerical results 
(black, dashed) for various values of the CP violating phase 5 at L — 2300 km. The average matter 
density for L — 2300km is p — 3.54 g/cm 3 . 



Ref. [31] to demonstrate the strength of their formula. The average Earth matter density 
for this baseline is 3.81g/cm 3 [35]. We consider the normal hierarchy case, Sm^ > 0, with 
the CP violating phase 5 set to zero. The differences between the exact and approximate 
formulae are plotted in Fig. 5(b). As can be seen, at this baseline, both the Asano-Minakata 
formula and our approximation work much better than the Cervera et al. formula. 

The comparison at a shorter baseline length of L = 810 km, which is the distance from 
Fermilab to NOi^A, is shown in Fig. 6. There, all three approximations work well, with 
our approximation being the most accurate. 

The situation changes at the longer baseline length of L = 8770 km, which is the 
distance from CERN to Kamioka [36], as can be seen in Fig. 7. There, the Cervera et 
al. formula greatly overestimates Piy^ — > v e ), while the Asano-Minakata formula leads to 
negative probability for E ~ 4 GeV. In comparison, our approximation remains accurate. 

The accuracy of our approximation for both the neutrino and anti- neutrino cases, and 
both mass hierarchies, for different values of the CP violating phase 5, is demonstrated in 
Figs. 8 and 9 for the two baselines L = 1300 km and L = 2300 km, respectively. These 



- 11 - 



distances correspond to those between Fermilab and Homestake (1300 km), and CERN 
and Pyhasalmi (2300 km) [37]. As is evident, our approximation maintains its accuracy 
for all energy ranges and mass densities. 



4 Applications 

4.1 Determination of the Mass Hierarchy from v e Oscillations 

Consider the v e survival probability in matter which is given by 

P{v e -)• v e ) 

= 1-4 |& e2 | 2 (l - |& e2 | 2 ) sin 2 \ - 4 |& e3 | 2 (l - |& e3 | 2 ) sin 2 \ 
+ 2|& e2 | 2 |[/ e3 | 2 ^4sin 2 ^sin 2 ^ +sinA 2 isinA 3 i^ 

i a II 12 Ia 12 12 \ ■ 2 ^21 . 2/ /)/ ^ • 2 ^31 

= l-4c 1 3S 12 (l-c 13 s 12 )sm — - sin (20 13 ) sin — 

, 12 ■ 2/ na l \ ( o • 2 ^21 . 2 ^31 , 1 ■ a A I 

+ s 12 sin (26> 13 ) I 2 sm -y- sin -y- + - sin A 2 i sin A31 I 

s 'i2~l n • 2/r>/)/ a/' ■ 2^21 . . 2^31 n . 2^21 . 2^31 1 ■ a ■ A ^ 

> 1 - sm (26 l 13 ) I sm — h sm — 2sm -ysin — sm A 2 i sm A31 I 

= 1- S in 2 (20' 13 ) sin 2 ^, (4.1) 

where we have assumed that a 3> Sm^i so that s' 12 ~ 1 is a good approximation. Similarly, 
we find: 



Piye 



A \tt \2\tt \2 • 2 ^21 . . i rr \2\tt 1 2 • 2 ^31 

4 \U e 2\ \U^ 2 \ sm z — + 4 |?7 e3 | z |C/ M3 | z sm z — 



+2 3? {u: z U^U e2 U; 2 ) ( 4 sin 2 ^ sin 2 ^ + sin A 21 sin A 31 



+4 J( eiM ) ( sin 2 — ^- sin A31 — sin 2 sin A21 ] (4-2) 



2 2 



A 12 12 I 12 2 x 12 12 2 r, / / / r\ • 2 ^21 ,4/2/22 -2 ^31 

4 si 2 ci 3 (c 12 c 23 + si 2 si 3 s 2 3 - 2si 2 c 12 Si 3 c 23 S23 cos J sm — + 4 s 13 c 13 s 23 sm — 
+2 s / 12 s / 13 c / 1 2 3S23 (c' 12 c 23 cos 5 - s' 12 s' 13 s 2 3) ^4 sin 2 sin 2 + sin A 2i sin A 3i ^ 

-4 s / 12 c / 12 s / 13 c / 1 2 3 s 23 c 23 sin 5 ^sin 2 sin A 3i - sin 2 ^y- sin A 2i ^ (4.3) 
> 4, sin 2 (2#' 13 ) sin 2 ^ 2 -, (4.4) 
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P(v e -)• V T ) 

= 4 \U e2 \ 2 \U r2 \ 2 sin 2 -|i + 4 |^ e3 | 2 |^ 3 | 2 sin 2 

+2 U (U* e3 U T3 U e2 U: 2 ) sin 2 sin 2 ^ + sin A 21 sin A 31 ^ 

sin — - sin A 3 i - sin — - sin A 2 i I 

a 12 12 I 12 2 , 12 12 2 o / / / c\ • 2 ^21 . /2 /2 2 -2 ^31 

= 4s 12 C 13 (c 12 S23 + s 12 S 13 c 23- 2s 12 c 12Sl3 s 23C23COS(5jsm — + 4 S 13 C 13 C 23 Sin — 

n I I 12 II x i ' ' \ I A "2 ^21 . 2 ^31 , • a • A I 

-2 s 12 s 13 c 13 C2 3 (c 12 s 23 cos <) + s 12 s 13 C2 3 j I 4 sm -ysm — — h sin A 2 i sin A 3 i I 

. , / / / /2 ■ r | ■ 2 ^21 . ~ . 2 A 3 i . ~ \ 

+4 Si 2 Ci2Si 3 Ci 3 s 23 C23 sm 6 I sm -y- sin A 3i - sin -y- sin A 2 i I 

s' 12 ~i 2 . 2 /not \ ■ 2 ^32 /, c n 
► c 23 sm (20 13 ) sm — . (4.5) 



From Fig. 3, it is clear that the factor sin 2 (2^ 3 ) in these expressions behaves quite differ- 
ently depending on the mass hierarchy. For normal hierarchy sin 2 (2^ 3 ) will peak around 
a ~ Sm 3 i but for the inverted hierarchy case it will not. This will become manifest if the 
factor sin 2 (A 3 2/2) also peaked at or near the same energy, 
For the normal hierarchy case, when a ~ <5rrt 31 we have 

5A 32 = X% - X'L ~ y/[X' + - (5m 2 31 + as 2 3 )} 2 + 4a 2 c 2 3 s 2 3 « 2s 13 a . (4.6) 

Therefore, 

A32 = 6Xv L ^ = 5 fp\ m 

2 4£ 2£ v ; Vg/ cm /V km / 
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(4.7) 



_ = 7T / pL 

2 V 54000 (km - g/cm 3 ) 

From Fig. 1, it is clear that pL < 54000 (km -g/cm 3 ) as long as the neutrino beam 
does not enter the core of the Earth, at which point the constant average matter density 
approximation breaks down. Therefore, in order to take A32/2 as close as possible to tt/2 
while preventing the beam from entering the Earth's core, we need L ~ 10000 km. 

For instance, if we take L = 10000km for which p = 4.53 g/cm 3 , we have pL ~ 
45300km • g/cm 3 . The value of A 32 /2 at resonance a ~ Sm^ is then 

vr 45300 , /Jn . 

- x = 0.42 vr , (4.8) 

2 54000 ' v ' 

leading to an oscillation peak/dip factor of sin 2 (A32/2) = 0.94. Using Eq. (2.2), the 
neutrino beam energy at which a ~ Sm^ is found to be 

E _ (Sml/eV 2 ) (2.47 x 10" 3 ) ^ ? _ 



GeV (7.63 x lO" 5 ) x (p/(g/cm 3 )) (7.63 x lO" 5 ) x (4.53) 

Indeed, in Fig. 10(a) we show the exact v e survival probabilities at L = 10000 km for 
both hierarchies, and we can see that the normal hierarchy case dips by over 95% around 
E ~ 6.5 GeV. Thus, our rough estimates give a highly reliable result. 

If we take a somewhat shorter baseline of L = 8770 km, which is the distance between 
CERN and Kamioka [36], we have p = 4.33g/cm 3 , and pL ~ 38000km • g/cm 3 . The value 
of A32/2 at resonance a ~ dm^ is then 

vr 38000 /t . 

- x = 0.35 vr , 4.10 

2 54000 ' v ' 

leading to an oscillation peak/dip factor of sin 2 (A32/2) = 0.8, which is still fairly prominent. 
Using Eq. (2.2), the neutrino beam energy at which a ~ is found to be 

E (JmjyeV 2 ) = (2.47 x 10' 3 ) ^ 

GeV " (7.63 x 10~ 5 ) x (p/(g/cm 3 )) (7.63 x 10" 5 ) x (4.33) ' ' 1 ' ' 

The actual oscillation peak occurs slightly off resonance around E = 6.5 GeV as can already 
be seen in Fig. 7. Comparison of P{v e — > u e ) at L = 8770 km with 5 = for the normal 
and inverted hierarchies are shown in Fig. 10(b). P{y e — > vS) is compared in Fig. 11(b). 

The differences between the normal and inverted hierarchies for both baselines is mani- 
fest, indicating that measuring these oscillation probabilities at this baseline would allow us 
to determine the mass hierarchy easily. (We consider the dependence on the CP violating 
phase 5 in the next section.) Eqs. (4.1), (4.4), and (4.5) also suggest that the measurement 
may provide a better determination of sin 2 623 ■ 

4.2 The "Magic" Baseline 

The "magic" baseline is the baseline at which the dependence of P{y e — > v^) on the CP 
violating phase 5 vanishes [38]. 4 Looking at Eq. (4.2), the only term without 5-dependence 



4 An illuminating discussion on the physical meaning of the "magic baseline" can be found in Ref. [39] 
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Figure 11. The dependence of the exact oscillation probability P(v e — > v^) on the CP violating 
phase 5 at (a) L = 7500km (p = 4.21g/cm 3 ), and (b) L = 8770km (p = 4.33 g/cm 3 ) for the normal 
(red) and inverted (green) mass hierarchies. 



is the |f/ e 3| 2 |^^3| 2 term. To make every other term vanish, we must have 



A 



sin 



21 



<5A 2 i T 
sin I L 



4E 



. 



Therefore, the magic baseline condition is 

5X 2 i 



\E 



-L 



nir 



n G 



(4.12) 



(4.13) 



If we are in the energy and mass-density range such that Sm^i < a < \Srn^i\, we can see 
from Fig. 4 that <5A2i ~ a = 2\/2GfN £ E, and the above condition reduces to 



V2G F N e L » 2nvr , 



(4.14) 



which is the usual magic baseline condition. Using Eq. (2.2), this condition for the n = 1 
case becomes 



^21 



that is 



L = (9.7 x 1(T 5 ) 



pL 



P 



g/cm 3 / \km J 



7T 



32000 . 



(4.15) 



(4.16) 



km • g/cm 3 

This is satisfied for L ~ 7500 km as can be read off of Fig. 1. Indeed, in Fig. 11(a) we plot 
the bands that P{y e — > Vu) at L = 7500 km sweeps for both mass hierarchies when 5 is 
varied throughout its range of [0, 2ir]. We can see that the dependence on 5 is very weak. 

However, if we look at Eq. (4.3) carefully, it is clear that all terms that include the CP 
violating phase 5 are multiplied by c' 12 which goes to zero when a 3> dm^i- Indeed, this 
was why 5 did not appear in Eq. (4.4). The condition a 3> Srn^i demands 



(- 



E 



g/cm 3 J \GeV 



> 1 



(4.17) 
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which is clearly satisfied around the oscillation peak for the L = 8770 km case just discussed 
in the previous section. Thus, P{y e — > v^) for this baseline is also only very weakly 
dependent on 5 as shown in Fig. 11(b). We can conclude that, in general, as long as 
Eq. (4.17) is satisfied, one does not need to be at a specific "magic" baseline to suppress 
the ^-dependence of P(y e — > u^). 

5 Summary 

We have presented a new and simple approximation for calculating the neutrino oscillation 
probabilities in matter. Our approximation was derived utilizing the Jacobi method [34], 
and we show in the appendix that at most two rotations are sufficient for approximate 
diagonalization of the effective Hamiltonian. The two rotation angles are absorbed into 
effective values of 0x2 an d #13. 

As explained in detail in the appendix, the approximation works when #13 = 0(e), 
where e = \/5rr^J\5m^\ = 0.17, a condition which has been shown to be satisfied by 
Daya Bay [23] and RENO [24]. Our formulae are compact and can easily be coded as 
well as be manipulated by hand for back-of-the-envelope calculations. The application of 
our method to finding the v e — > tv, v T resonance conversion condition, and that to the 
determination of the 'magic' baseline [38, 39] have been demonstrated. 

In this paper, only the matter effect due to Standard Model W exchange between 
the neutrinos and matter was considered. New Physics effects which distinguish between 
neutrino flavor would add extra terms to the effective Hamiltonian, which would require 
further rotations for diagonalization. This has been discussed previously in Ref. [33], and 
the potential constraints on New Physics from long baseline neutrino oscillations experi- 
ments in Refs. [40-42]. Updates to these works will be presented in future publications 
[43]. 
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A Conventions, Notation, and Basic Formulae 



Here, we collect the basic formulae associated with neutrino oscillation in order to fix our 
notation and conventions. 



A.l The PMNS Matrix 

Assuming three-generation neutrino mixing, the flavor eigenstates \v a ) {a = e,fj,,r) are 
related to the three mass eigenstates \uj) (j = 1,2,3) via the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) matrix [44-46] 



(^PMNs)aj = (Va\Vj) , 



that is, 



\v a ) {v a \vj) 



Wat) = ^2 ^ ( u j\ u a) 
3=1,2,3 

The standard parametrization is given by 



(VpMNS)aj Wa) , 

Q=e,fi,r 
3=1,2,3 



Wl 



th 5 



VpMNS = UV , 



(A.l) 



(A.2) 



(A.3) 



u = ^23(^23,0) Ri 3 (e 13 ,6) R 12 (0 12 ,0) 



V 



1 

C 2 3 «23 
-S23 C 2 3 

C12C13 



1 

-si 3 e lS ci 3 

S12C13 

AS 



Cl2 3X2 
-SX2 C12 

1 



«i3e" 



i5 



S12C23 - Cl2Sl3«23e C12C23 - Si 2 Sl3S23e C13S23 



Sl2«23 - Cl2Sl 3 C 2 3e 

diag(l,e 



id 



-C12S23 - Sl2Sl 3 C23e J<5 C13C23 



Ja 2 i/2 g «a 3 i/2\ 



(A.4) 



Here, Rij(9, 5) denotes a rotation matrix in the ij-plane of clockwise rotation angle 9 with 
phases ±<5 on the off-diagonal ji and ij-elements, respectively, and Sij = sin9ij, Cij = 
cos 9{j. Without loss of generality, we can adopt the convention < 9{j < 7r/2, < 6 < 2tx 
[47]. Of the six parameters in this expression and the three neutrino masses, which add up 
to a total of nine parameters, neutrino— meutrino oscillations are only sensitive to six: 

• the three mixing angles: 9yi, 023, #13, 

• two mass-squared differences: Srri^x, Sm^, where 5mfj = mf — m|, and 

• the CP-violating phase: 5. 

5 Cervera et al. in Ref. [22] use a different convention in which the sign of 8 is flipped. 
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The Majorana phases, CK21 an d «3i, only appear in lepton- number violating processes 
such as neutrinoless double beta decay, and cannot be determined via neutrino— ^neutrino 
oscillations. The absolute scale of the neutrino masses also remain undetermined since 
neutrino oscillation is an interference effect. 

A. 2 Neutrino Oscillation 

If a neutrino of flavor a is created at x = with energy E, then the state of the neutrino 
at x = is 



Wa, {x = 0)) = K) = ^(VpMNs)aj Wj) ■ 

3=1 



At x = L, the same state is 
3 

\v afl (x = L)) = (^MNS ) aj n 

3=1 

Assing uij CBwe can approximate 



(A.5) 



• a .\ Vj ) = e^^e^-^ypMNs)^) • (A.6) 
3=1 



so that 



and we find 



Pj 



Pj ~ Pi 



E 2 -m 2 = E 



' 2E ' 



2E 



+ 



Sin?! = m 2 — m\ , 
K,o(x = L)) = e^J^exp -i-^L (V PMNS )* aj \ Vj ) 



3=1 



(A.7) 
(A.8) 

(A.9) 



Therefore, the amplitude of observing the neutrino of flavor /3 at x = L is given by (dropping 
the irrelevant overall phase) 



Ap a = {vp\v afi (x = L)) 



("hi (VpMNs)/3k 

.k=l 

3 / 



^exp| 

■3=1 
Sm 2 ^ 



5m 2 n \ 
' i ^^ L ) (^PMNs)aj Wj) 



= ^(^PMNs)/3i exp( -i^-L ] (VpuNS) aj 



3=1 
3 



Pj exp -1 



3=1 \ 
~U exp(-i^L ) I/* 



2E 



Pa 



(A.10) 
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where 



5M' 




5m\ x 
5m 2 31 



and 



H n = U 5M 2 U ] . 



(A.ll) 



(A.12) 



Thus, the probability of oscillation from \v a ) to \vp) with neutrino energy E and baseline 
L is given by 



P(v a ->■ z^g) 




2E 



A 



4 ^ R {KiUpiU^) sin 2 -^ + 2^9 [U*iU pi U aj U^) sin A 



(A.13) 



where 6 
A 

Since 



5m 2 , 

—L 

2E 



2.534 



5ml 



L 



eV 2 J\ E J \km 



5m: 



m,-, 



A 



32 



A 



31 



A 



21 



(A.14) 



(A.15) 



only two of the three Ay's in Eq. (A.13) are independent. Eliminating A32 from Eq. (A.13) 
for the a = j3 case yields 



P{v a ^v a ) = l-A\U a2 \ 2 {l-\U< 



+ 2 |£/a2| 2 |£^a3| 2 



a2\ ) Sin 

A 



^21 



Mu a3 \ 2 {1 - \U a3 \ 2 ) 



2\ -2 ^31 



sm 



4 sin 2 



21 



A 



sm 



31 



+ sin A 2 i sin A31 



(A.16) 



and for the a 7^ (5 case we have 



A 



P{y a -»• vp) = MU a2 \ 2 \Up 2 \ Z sin 2 -^ + 4|C/ Q3 | 2 |^ 3 | 2 sin 



A 



31 



+2 SR {U* 3 Up 3 U a2 U; 2 ) (^4 sin 2 ^ sin 2 ^ + sin A 21 sin A31) 



+4 J, 



A 



sm 



21 



sin A31 — sin' 



A 



31 



sin A 



21 



(A.17) 



where J( a B) is the Jarlskog invariant [49]: 



+Z(U* 1 U fn U a2 U, 



B2i 



+Z(U* 2 Up 2 U a3 U 



83) 



+%{K 3 u f}3 u al u f 



BlJ 



6 Note that our notation differs from that of Cervera et al. in Ref. [22]. There, the symbol Aij is defined 
without the factor of L, that is, Aij — 5m1j/2E. It also differs from that used by Freund in Ref. [17] where 
A = 8m% 1 , and A = Sm^iL/4E. Huber and Winter in Ref. [38] define A = Sm^L/AE, which is also used 
in Ref. [48]. So care is necessary when comparing formulae. 
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133) 



MK 3 u 03 u a2 u, 



132/ 



In the parametrization given in Eq. (A. 4), we have 



J 



(e,M) 



J 



(e,r) 



(r,e) 



J (r,M) 



J, 



JsinJ 



with 



J 



s 12 c 12 s 13 c 13 s 23 c 23 



(A.18) 

(A.19) 
(A.20) 



The oscillation probabilities for the anti-neutrinos are obtained by replacing U a i with its 
complex conjugate, which only amounts to flipping the sign of 5 in the parametrization of 



Eq. (A.4). It is clear from Eq. (A. 16) that P(u c 



which is to be 



expected from the CPT theorem. For flavor changing oscillations, only the Jarskog term 
in Eq. (A. 17) changes sign. 

A.3 Matter Effects 

If the matter density along the baseline is constant, matter effects on neutrino oscillations 
can be taken into account by replacing the PMNS matrix elements and mass-squared 
differences with their "effective" values in matter: 



Un 



Ur 



where U is the unitary matrix that diagonalizes the modified Hamiltonian, 



(A.21) 



H„ 



U 



Ai 





" 


~t 

u = u 


"o 




"a 0" 










5m 2 2l 


17+ + 











A 3 




Smli 








(A.22) 



5M 2 



Hn 



and 



A 



5X 



2E 



L 



'.i 



Ai — A, 



(A.23) 



The factor a is due to the interaction of the \u e ) component of the neutrinos with the 
electrons in matter via ^-exchange: 



a = 2V2 G F N e E 



(A.24) 



Assuming N e = N p ~ N n in Earth matter, N e for mass density per unit volume of p can 
be expressed using Avogadro's number Na = 6.02214129 x 10 23 mol _1 as 

^23 /„ TV ,3\ I P 



pN A /2 = (3.011 x W Z6 /cm 6 ) x 



(A.25) 



v g/cm 3 

Thus, putting back powers of he to convert from natural to conventional units, we find 

a = 2V2 G F N e E x (he) 3 
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(7.63 x l(T 5 eV 2 



(- 



E 



(A.26) 



v g/cm 3 y \^GeV / 

For anti- neutrino beams, a is replaced by —a in Eq. (A. 22). Note that a is -E-dependent, 
which means that both U and A™ are also independent. It is also assumed that E <C Myy 
since the VF-exchange interaction is approximated by a point-like four-fermion interaction 
in deriving this expression. 



B Jacobi Method 



B.l Setup 



As mentioned in the introduction, it is possible to write down exact analytical expressions 
for A{j and U a i [12]. However, simpler and more transparent approximate expressions can 
be obtained using the Jacobi method as will be shown in the following. 
We introduce the matrix 

Q = diag(l, l,e iS ) , (B.l) 
and start with the partially diagonalized Hamiltonian: 

H' a = Q^H a UQ 


8m 2 2l + U j 


8m 2 2l Q + aQ 





Q 



8ml x 



U \ Q 







8ml x 





5m 2 21 







+ a 



2 „2 



a 



U* el U el U* el U e2 

u* e2 u el u* e2 u e2 u: 2 u e3 

U* eS U el U* e3 U e2 U* e3 U e3 



Q 



Ci oC 



12 L -13 



Cl2Sl2C 13 C12C13S13 

C12S12C13 S12C13 S12C13S13 



C12C13S13 S12C13S13 
2 



'13 



ac 12 c 13 OCi 2 Sl2Cl 3 aCi 2 Ci 3 Si3 

aci 2 s 12 cj 3 as\ 2 c\ 3 + bm\\ asi 2 c 13 si 3 



aci 2 ci 3 si3 asi 2 ci 3 si3 asf 3 + Sm^ 



(B.2) 



The matrix Q serves to rid H' a of any reference to the CP violating phase 5. The strategy 
we used in our previous papers [32, 33] was to approximately diagonalize H' a through the 
Jacobi method using 



5m 2 



21 



Mil 



0.17 



(B.3) 



as the parameter to keep track of the sizes of the off-diagonal elements. We argued that 
approximate diagonalization was achieved when the off-diagonal elements were of order 
0(e 2 s 13 \8m 2 l \). 

Note that our e differs from Asano and Minakata's e in Ref. [31] where 



5m 2 



21 



Mil 



0.03 . 



(B.4) 



That is, e = e 2 . So care is necessary when comparing formulae. 
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B.2 Diagonalization of a 2 x 2 Matrix 



Recall that for 2x2 real symmetric matrices, such as 



M 



a f3 
P 7 



diagonalization is trivial. Just define 



R 



and we obtain 



with 



where 



COSCJ 



Ai 

A-2 



R^MR 



a, /3,j e 



(B.5) 



tan 2uj 



7 — a 



Ai 
A 2 



acr. 



7 4 _ _ ( a + 7) T ^(a - l) 2 + 4/? 2 



7C U 



as;. 



(a + 7) ± - 7 ) 2 + 4^ 2 



, (B.6) 



(B.7) 



(B. 



where the upper and lower signs are for the cases a < 7 and a > 7, respectively. The 
Jacobi method [34] entails iteratively diagonalizing 2x2 submatrices of a larger matrix 
in the order that requires the largest rotation angle at each step. In the limit of infinite 
iterations of this procedure, the matrix will converge to a diagonal matrix. 

In the case of H' a given in Eq. (B.2), at most two iterations are sufficient to achieve 
approximate diagonalization, neglecting off-diagonal elements of order 0(s 2 s^ldm 2 ^) , re- 
gardless of the size of a. We demonstrate this in this appendix. 



B.3 Neutrino Case 

B.3.1 Mixing Angles and Mass-squared Differences 

Let us first evaluate the sizes of the sines and cosines of the three vacuum mixing angles 
in comparison to the parameter e defined in Eq. (B.3). The current best fit values for the 
mass-squared differences and mixing angles are listed in Table 1. The sines and cosines of 
the central values of the mixing angles are 

S23 = 0.71 , C 23 = 0.71 , 

512 = 0.55 , C12 = 0.84 , (B.9) 

513 = 0.15 , C13 = 0.99 . 

Therefore, S13 is O(e) while all other sines and cosines are 0{\). 
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Figure 12. (a) The dependence of 6' 12 on j3 — — log e (a/j^m^l). (b) The /3-dependence of A±. 



B.3.2 First rotation 

The effective hamiltonian we need to diagonalize is 



ac 12 c 13 ac 12 si 2 cf 3 aci 2 ci 3 si 3 
ac 12 s 12 cj 3 as\ 2 c\ 3 + bra\ x as 12 c 13 s 13 
ac 12 c 13 s 13 asi 2 ci 3 si3 asj 3 + 5m 2 3l 
aO(l) aO(l) aO{e) 
aO(l) aO(l) + 5ml x aO{e) 
aO{e) aO{e) aO{e 2 ) + 5m 2 31 



(B.10) 



Of the off-diagonal elements, the 1-2 element is the largest regardless of the size of a. 
Therefore, our first step is to diagonalize the 1-2 submatrix. 








y"^ a sin[# 12 '] 






a cos[#i 2 '] 



13 



(a) 



(b) 



Figure 13. (a) The dependence of s' 12 = sm9' 12 and c' 12 = cos6*' 12 on (3 — — log e (a,/ 1 <5tti.§ 1 |) . (b) 
The dependence of as' 12 and ac' 12 on /?. The values are given in units of l^mlJ. The asymptotic 
value of ac[ 2 is Sm^s^c^/c^ « 0.014 l^m^j = O^^m^]). 
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Define 



V 



s<p 



1 



(B.ll) 



where 



COS (f , Su 



sin ip , tan 2ip 



ac\ 3 sin 2$ 



12 



c> m 2i — ac 2 3 cos 26>i2 



< (f < 



7T 



(B.12) 



Using V, we find 

H" = V^H'V 











aC 12 Cl3Sl3 
as' 12 Cl3Sl3 



ac' 12 ci 3 si3 as' 12 ci3Si3 as? 3 + tfraj; 



where 



and 



c 'l2 



cos ( 



12 ' 



i 

s 12 



sin ( 



12 ) 



? 12 = 0l2 + f 



{5m 2 2l + acf 3 ; 



21 



ac 



f 3 ) 2 + Aac 2 13 s 2 2 5m 2 



21 



* 2 
The angle 9' 12 = #12 + </> can be calculated directly without calculating ip via 

5m 21 sin 2^i2 



tan 20' 12 



o~m 21 cos 2^12 — ac 2 3 



^6*12 < 0'i2 < 



(B.13) 

(B.14) 
(B.15) 

(B.16) 



The dependences of 6' 12 and A'_j_ on j3 = — log e (a/\Sm 2 i\) are plotted in Fig. 12. Note 
that 0' 12 increases monotonically from 9\ 2 to ir/2 with increasing a. The /3-dependence of 
s' 12 = sin 6^2 and c' 12 = cos6" 12 are shown in Fig. 13(a). For a ^> 5m 21 , s' 12 and c' 12 behave 
as 



s 'l2 



c 'l2 



1 



s 2 c 2 



12°12 



5m 2 21 



ac 



13 



+ 



•S12C12 



5m 2 2l 



ac 



13 



+ Sl2Ci 2 (c 



12 



4) 



ac 



13 



+ • 



(B.17) 



Therefore, for a 3> o~m 21 we have as' 12 « a while ac' l2 rzs (5m 21 si2Ci2/c 2 3 = e 2 \5m 2 1 \si2Ci2/c 2 3 
0.014 |<5m 31 | = 0(e 2 \5m 2 1 \). This behavior is shown in Fig. 13(b). Note that ac' 12 never 
grows larger than 0(e 2 |5m 31 |) for any a. 

The values of A^_ away from the level crossing point a ~ 5m 21 for the a <C Sm 21 case 
are given by 



\/ _ 22 
A_ — ac 13 c 12 



1 



s 12 



ac 



13 



5m 2 21 



'12 



(C?2 



-fa) 



13 



5m 2 21 



+ ••■ 



A' 



(5m 21 



1 + 4 



13 



^21 



,22 

"T s 12 c 12 



or 



13 



^21 



+ s 12 c 12( c 12 ~~ s 12) 



ac 



13 



5m 2 2l 



+ 



(B.18) 



- 24 - 




Figure 14. (a) The /3-dependence of <f> for the normal and inverted hierarchies, (b) The /in- 
dependence of the difference 0' — <p. 



and those for the a 3> 5m 21 case by 



A' 



A' 



s 2 
s 12 



(57712! 



Of 



13 



s 12( c 12 ~~ s 12) 



^21 



acr 



13 



+ 



1 + s 



12 



5m 2 2l 



ac 



1.3 



,22 
' s 12 c 12 



5m 2 2l 



ac 



1.3 



,2 2 / 2 2 \ 

' s 12 c m c 12 ~~ S 12J 



5m 2 2l 



ac 



13 



+ 



B.19) 



We will use this expansion for A+ later. Thus, the asymptotic values of A± are A'_ 



nr 2 r 2 \' 
Mc 13 c 12> A + 



5m 2l in the a — > limit, and A'_ 



5m 



21 L 12 



ac 2 3 in the a — > oo limit. 



B.3.3 Second rotation 

The effective hamiltonian after the first rotation was given by Eq. (B.13). When a < 5m 21 , 
both non-zero off-diagonal elements are of order 0(ea) < 0(e 3 \5m 2 1 \), since s' 12 and c' 12 
are both O(l) in that range as can be discerned from Fig. 13(a). However, as a increases 
beyond Sm 2 2 an d 0' 12 approaches 7r/2, we have as' 12 — > a, ac' 12 — > 0{e 2 \5m\ l \) and the 2-3 
element becomes the larger of the two. Therefore, a 2-3 rotation is needed next. 
We define 



W 



where 



C0 



cos < 



S<j> 



sin ( 





tan 2(f) 



as' 12 sin 26>i3 



^ m 3i + as 



13 



Ai 



(B.20) 



(B.21) 



The angle 

l 31 



is in the first quadrant when 5m 2 1 > 0, and in the fourth quadrant when 



5m 2 i < 0. Then, 



H'" = W ] H"W 



-aC^CisSisS^ 
ac' 12 Cl3Sl3C0 



-ac' 12 ci3Si 3 S0 ac'^cissisc^ 





4 



(B.22) 
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where the upper(lower) sign corresponds to the normal(inverted) hierarchy case with 



[ A'+ + (6ml, + as 2 13 ) ] ± J [ X' + - (Sm 2 31 + o«? 3 ) ] + 4a 2 



b 12 L 13 *13 



,// _ L T V OX W/ J V L -T V OJ. 10/ J X4 xo xo ,_. , 

A ± = . [15.26) 

The /3-dependences of A± and 4> are shown in Fig. 4 ((a) and (b)), and Fig. 14(a), respec- 
tively, for both mass hierarchies. For the normal hierarchy case, Sm 31 > 0, the values of 
\± away from the level crossing point a ~ Sm^i are approximately 

A" « 5ml, + as i3 i 

A'l A'+ , (B.24) 



when a <C (fon^i, and 



A+ W O + S?3 5m 31 + c 13 s 12^ m 21 > 
2 r_2 1 „2 „2 r_2 



A" « c 13 <5m 31 + s 13 s 12 (5m 21 , (B.25) 

when a 3> 8m\,. For the inverted hierarchy case, <5m 31 < 0, where there is no level crossing, 
the values of A" are approximately 

A" rj Jmli < , A" « A+ , (B.26) 

for all a. 

At this point, we argue that the angle <p defined in Eq. (B.21) is well approximated by 
the angle <f>' which we define via 

t 2(f)' ~ a sin 2^13 (B 27) 

(5m 2 -, — 5m 2 , 1 sf 2 ) — a cos 26 , i 3 

This approximation is obtained by first noting that <fi is significantly different from zero 
only when a S> 5m 2 ,. The expansion of X', in the denominator of the right-hand-side of 
Eq. (B.21) in powers of <5m 21 /a was given in Eq. (B.19). Keeping only the first two terms, 
and noting also that s' 12 ~ 1 to the same order when a 3> 5m|i (c.f. Eq. (??)) we obtain 
Eq. (B.27). The /3-dependence of the difference </J — <j> is plotted in Fig. 14(b), and we can 
see that the disagreement is at most 0(e 4 ). Thus, we replace (j) with <\>' in the following. 

Now, the effective Hamiltonian after the second rotation was given by Eq. (B.22). Note 
that all of the non-zero off-diagonal elements include the factor ac' 12 , which is never larger 
than 0(e 2 |<5r7T,| 1 1 ) regardless of the value of a as discussed above. They also all include 
a factor of si 3 , which is 0(e) as we have seen in Eq. (B.9). Therefore, all off-diagonal 
elements of H"' are of order 0(e 2 si^\5m'l 1 \) = 0(e 3 \5m 2 1 \) or smaller regardless of the 
size of a. Note that had the value of si 3 been smaller, the sizes of the neglected terms 
would have been proportionately smaller also. We conclude that, at this point, off-diagonal 
elements are negligible and a third rotation is not necessary. 
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B.3.4 Absorption of <fi' into #13 

From the above consideration, we conclude that the matrix which diagonalizes H' a , Eq. (B.2), 
is given approximately by VW, and that the effective neutrino mixing matrix becomes 



U « UQVW = R2 3 (9 23 ,Q)R 13 (9i 3 J)R 12 (9 12 ,0)QRi 2 (v,0)R 23 (<t>\Q) 

V 



(B.28) 



W 



Using 



we find 



U 



R 12 (6 12 ,0)Q = QRi 2 (9 12 ,0) , 
Ri3(0 13 ,5)Q = QR 13 (9 13 ,Q) , 



U « R2 3 (e 23 ,0)Ri 3 (0 l3 ,5)R 12 (0 l2 , 0)QRi 2 {<p, 0)R 23 {<f>' ,0) 
= #23(^3, 0)QR 13 (9 13 , 0)Ru(0j3, d)Ri 2 (if, 0)^23(0', o) 
= #23(^23, 0)QR 13 (9 13 , 0)R 12 (e 12 + if, 0)Rw(i>', 0) 

= R23(023,0)QRl3(0 13 ,0)R 1 2{e' 12 ,0)R 23 ((j ) ',0) . 



(B.29) 



Here, we argue that 



i2i2(0' 12 , 0)^23(^,0) « Ri 3 (cf>',0)Ri 2 (e' l2 ,0) 



(B.30) 



(B.31) 



that is, the 2-3 rotation becomes a 1-3 rotation when commuted through Ri 2 (9' 12 ,0). This 
is due to the fact that 4>' only becomes non-negligible when a 3> 5m^ 2 where s' 12 ~ 1 and 



'12 



0, which means 



Ri2(e' 12 ,o) 

and it is straightforward to see that 



1 
-10 
1 



(B.32) 
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(B.33) 


1 












1 





where = sin<// and = coscp'. In the range o < 5m 21 , the angle is very small and 
both .R23 (</>', 0) and i?i3(0',O) are approximately unit matrices and Eq. (B.31) is trivially 
satisfied. Curiously, this approximation breaks down around a ~ for the normal 

hierarchy case when #13 is 0(e 2 ) or smaller, as is discussed in appendix C. However, given 
that the current experimentally preferred value of #13 is O(e), the approximation is valid. 
Thus, 

U « #23(023, 0)2^13(013, 0)Rl2(9' 12 ,Q)R2 3 (<j>',0) 

« #23(^23, 0)QR 13 (9 13 , 0)R 13 (<p', 0)R 12 (9' 12 , 0) 
= i? 23 (023, 0)Q i?i 3 (0i3 + <t>', 0)R 12 (9' 12 ,0) 
= R23(9 23 , 0)QR 13 (9' 13 , 0)R 12 (9[ 2 , 0) 
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R23(023,o)R 13 (e' 13 ,s)R 12 (e' 12 ,o)Q, 



where we have defined 



#'l3 = 013 + <t>' ■ 

This angle can be calculated directly without calculating 4>' via 

((knjj — 5m2isf 2 ) sin 2^i3 



tan 26>i 3 



(dm^ — dm^s^) cos 2^i3 — a 



(B.34) 
(B.35) 

(B.36) 



The diagonal phase matrix Q appearing rightmost in the above matrix product can be 
absorbed into the redefinition of the major phases and can be dropped. Thus, we arrive at 
our final approximation in which the vacuum mixing angles are replaced by their effective 
values in matter 



Cl2 
#13 

5 



#12 
#13 
#23 

5, 



= #12 + <P , 
= #13 + , 



(B.37) 



and the eigenvalues of the effective Hamiltonian are given by 



Ai 

A 2 
A 3 



A'_, 



(B.38) 



Note that of the mixing angles, only #12 and #13 are shifted. #23 and S stay at their vacuum 
values. 

B.4 Anti-Neutrino Case 
B.4.1 First Rotation 

For the anti-neutrino case, the matter effect parameter a acquires a minus sign. Thus, the 
effective hamiltonian to be diagonalized is 



2 2 
— ac 12 c 13 



aci 2 si2q 3 -aci 2 ci 3 si3 

aCi2Si 2 cf 3 -aS?2 C 13 + ^ m 21 _as 12Cl3Sl3 

2 



^ m 31 



-aci2Ci3Si3 -asi2Ci3Si3 -as 13 
-aO{\) -aO(l) -aO(e) 
-aO(l) -aO(l) + 5m 2 2l -aO(e) 
-aO(e) -aO(e) -aO(e 2 ) + Sm 2 31 



(B.39) 



The largest off-diagonal element is the 1-2 element. Therefore, our first step is to diagonalize 
the 1-2 submatrix. 



Define 



V = 



C(p S(p 

-s v c v 
1 



(B.40) 
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(a) 



(b) 



Figure 15. (a) The dependence of 9 12 on j3 — — log e (a/|5m| 1 1). (b) The /3-dependence of A ± . 



where 



COS if 



j °<p 



sm<p , t&n2(p 



acf 3 sin 26?i2 



5m 21 + acf 3 cos 2#i2 



7T 



< 99 < 



Using V we find 



h" = v ] h'v 



A_ U -ac 12 ci 3 si 3 

A + -0813013513 

-ac' 12 ci 3 si3 -as' 12 ci3Si3 -asf 3 + Jm^ 



where 



and 



'12 



A ± = 



cos ( 



12, s 12 = Sln #12, #12 = 12 + 

(5m|i - ac 2 3 ) ± V (5m 2 i + ac^ 3 ) 2 - Aac{ z s'{ 2 5m 2 2l 



The angle # 12 can be calculated directly without going through ip via 

5m 21 sin 2#i2 



tan2# 12 



5m 21 cos 26*12 + oc? 3 ' 



< 



12 ^ $12 



(B.41) 

(B.42) 

(B.43) 
(B.44) 

(B.45) 



The /3-dependences of 12 and A_|_ are shown in Fig. 15. Note that in contrast to the 
neutrino case, there is no level crossing. 9 l2 decreases monotonically toward zero as a is 
increased. For a S> Sm 2 i, s' 12 and c' 12 behave as 



S\2 = S12C12 



5m 2 



21 



a.c 



13 



Sl2Cl 2 (cf 2 - s\ 2 ) 



5m\ x x 1 



nc 



c 12 



s 2 2 
6 12 c 12 



5m 



2 \2 
21 



13 



+ • 



a a 



13 



+ 



(B.46) 



and we see that, this time, we have ac' 12 « a and as' l2 « 5m 21 si2Ci2/ 'cf 3 = 0(e 2 |(5m 31 | 
These /3-dependences of s' 12 , c' 12 , as' 12 , and ac' 12 are shown in Fig. 16(a) and (b). 




Figure 16. (a) The /3-dependence of s' 12 = sin# 12 and c' 12 = cos6*' 12 . (b) The /3-dependence of as' x 
and ac' 12 . The asymptotic value of as' 12 is <5m 21 Si 2 Ci 2 / 'cf 3 0.014 l^m^l = O (e 2 1 <5ttt,§ 1 | ) . 



In the range a <C Sm^, A_|_ can be expanded as 



T 7 2 2 

a_ — — ac 13 c 12 



2 / ac 

l + sj 2 



ac 



13 



(<i('j;j \ 2/2 2 \ ( 



+ 



5m 2 



21 



'12 



ac 



13 



2 „2 



+ s 12 c 



12^12 



1.3 



2 2 / z 
s 12 c m c 12 ~~ s 



L) 



ac 



13 



+ 
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while in the range a 3> Sm^, we obtain 



-ac 



13 



<5mJi c i2 



'12 



ac 



.22 
+ • s 12 c 12 



1 + S 



12 



1.3 

5m|i 



<5™21 



ac 



13 



ac 13 

s 12( c 12 ~~ s 12 



2 2 
s 12 c 12 



^21 



ac 



13 



H 



ac 



13 



+ 



The asymptotic values are thus A_ 
A + c^Smlx, A_ 



■ Sm^, A^ 
-acf 3 , in the limit a — > oo. 



-acf 3 cf 2 , m the limit a 



(B.48) 
0, and 



B.4.2 Second Rotation 

After the first rotation, the effective hamiltonian was given by Eq. (B.42). When a < Sm^i, 
both non-zero off-diagonal elements are of order 0(ea) < 0(e 3 |<5m 31 |). In contrast to the 
neutrino case, as a increases beyond £m 21 , the angle 12 approaches 0, and it is the 1-3 
element that becomes the larger of the two. Therefore, a 1-3 rotation is needed next. 
We define 

C0 

IT - 10 , (B.49) 

-Scf, C0 
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where 



c<b 



COS ' 



COS ( 



tan 20 



ac' l2 sin 2#i3 



5m 



(B.50) 



31 



as 



13 



The angle is in the fourth quadrant when 5m 31 > 0, and the first quadrant when 5m 31 < 
0. Using W, we find 



h"' = w ] h"w 







as 12 ci 3 si 3 S0 
as 12 ci3Si3S0 A + -as 12 ci3Si3C</, 

-as' 12 ci3S 13 C0 A ± 



(B.51) 



where the upper(lower) sign corresponds to normal(inverted) mass hierarchy with 



[ A'_ + (^mli - oa? 3 ) 



A_ — (5m 



as 



13J 



(B.52) 



The /3-dependence of A ± and 4> are shown in Fig. 4 ((c) and (d)), and Fig. 17(a), respectively, 
for both normal and inverted mass hierarchies. For the normal hierarchy case, 5rn 3l > 0, 
there is no level crossing, and A_j_ are well approximated by 



A, 



A" 



A' 



(B.53) 



Level crossing occurs for the inverted hierarchy case, 5m 31 < 0, in which we have 



-5m 



31 



a»13 > 



(B.54) 



when a <C <5"2-3i, and 



A 



-c 2 13 5m 2 31 + 4si2^2i , 
-a - sfg^mli + Ci 3 si 2 5m2i , 



(B.55) 



when a S> ^mf^. 

Here, as in the neutrino case, we approximate eft with the angle 



defined via 



tan 20 



a sin 26 



13 



(5m? 



31 



sf 2 5m 21 ) + a cos 2 2#i3 



(B.56) 



which is obtained by using Eqs. (B.46) and (B.48) on Eq. (B.50). The difference between 
4> and 4> is shown in Fig. 17(b), and it is clear that the difference is negligible. 

Now, the effective Hamiltonian after the second rotation was given by Eq. (B.51). Note 
that all of the non-zero off-diagonal elements include the factor as' 12 , which is never larger 
than O^l^m^D regardless of the value of a as discussed above. They also all include 
a factor of si 3 , which is 0(e) as we have seen in Eq. (B.9). Therefore, all off-diagonal 
elements of H a are of order 0(e 2 si 3 \5m 2 l \) = 0(e 3 |<5r?i3 1 |) or smaller regardless of the 
size of a. We conclude that, at this point, off-diagonal elements are negligible and a third 
rotation is not necessary. 
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Figure 17. (a) The /3-dependence of </> for the normal and inverted hierarchies, (b) The f3- 
dependence of the difference — <p. 



B.4.3 Absorption of <fi into #13 

From the above consideration, we conclude that the matrix which diagonalizes H' a , Eq. (B.39), 
is given approximately by VW, and that the effective anti-neutrino mixing matrix becomes 



U « UQVW 



#23(023, 0)i2 13 (0i3, S)Ri 2 (e 12 , Q\QR 12 (<p, 0)RM , 0) 

V 



(B.57) 



U 



w 



As in the neutrino case, we find 



U « R23(023 , 0)R 13 (e 13 ,5)R 12 (e 1 2,0)QRi2&,0)R 13 (^,0) 

= #23(023, o)QR 13 (e 13 , o)R 12 {e 12 , o)r 12 &, o)r 13 {4>', o) 

= #23(023, O)Q#1 3 (013, 0)^12(012 + <f>, 0)R 13 $ ', 0) 
= #23(023, 0)QR 13 (9 13 , 0)R 12 {rf 12 , 0)R 13 $,0) ■ 



Here, we argue that 



#12(^12, 0)#13 (0,0) « #13(0,0)^12(012,0), 



that is, the 1-3 rotation passes through i?i2(0i2,O). This is due to the fact that 
becomes non-negligible when a 3> Sm1 2 where s' 12 ~ and c' 12 ~ 1, which means 



#12(012, 0) 



1 
1 
1 



(B.58) 



(B.59) 
ft only 



(B.60) 



thus any matrix will commute with i?i2(0i2, 0). In the range a < Sm^i, the angle 4>' is very 
small and both R 23 ((j>', 0) and Ri 3 (4>',0) are approximately unit matrices and Eq. (B.59) is 
trivially satisfied. The accuracy of this approximation is discussed in appendix C. There- 
fore, 



U « #23(023,0)2 #13(013, O)#1 2 (012,O)#1 3 (0,O) 
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#23(023, 0)QR 13 {6 13 ,0)R l3 {cf> , 0)R l2 (e i2 ,0) 

#23(023, 0)2^13(013 + 4>', 0)Rl 2 {0' 12 , 0) 
#23(^23, 0)Q #13(^3, 0)Rl2(e' 12 , 0) 

#23(023, 0)Ri30' 13 ,S)R 12 (9[ 2 ,O)Q , (B.61) 



where we have defined 

013 



013 + 4>' ■ (B.62) 



This angle can be calculated directly without calculation <f>' via 

tan 24 = jY 31 ~/f 2 f' 2)Sin J 613 ■ CB.63) 
(ira^ — om2 X s\ 2 ) cos 2^13 + a 

The phase matrix Q appearing rightmost in the above matrix product can be absorbed 
into the redefinition of the major phases and can be dropped. Thus, we arrive at our final 
approximation in which the vacuum mixing angles are replaced by their effective values in 
matter 

?12 + V , 
+ , 

(B.64) 

and the eigenvalues of the effective Hamiltonian are given by 
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+ 012 = 






013 ~ 


* 013 = 


023 ~ 


* 023 , 
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Ai f 




\ 2 f 




A 3 f 





(B.65) 

Note that of the mixing angles, only 9± 2 and #13 are shifted. 9 23 and 5 stay at their vacuum 
values. 



C Commutation of i?i3 and R 23 through R\ 2 

In the derivation of our approximation formulae above, Eqs. (B.31) and (B.59) played 
crucial roles in allowing the second rotation angle to be absorbed into #13. In this appendix, 
we evaluate the validity of these approximations. 

C.l Neutrino Case 

The difference between the two sides of Eq. (B.31) is given by 



5R = R 12 (9' l2 ,0)R 23 (^,0)-R l3 ( ( f>',0)Ri 2 (9[ 2 ,0) 
^12(1-0 _(i_ s ' 12 )^' 







C 12 s <f> 



-C 
-(1 



12 V 1 -C^) 



C 12 s <f> 



'12^0 







(C.l) 
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(b) Inverted Hierarchy 



Figure 18. /3-dependence of c' 12 , 1 — s' 12 , sL and 1 — for (a) normal and (b) inverted hierarchies. 
The behaviors of c' 12 and 1 — s' 12 are common to both. 
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Figure 19. /3-dependence of the non-zero elements of 5R for (a) normal and (b) inverted hierarchies. 



It is clear that 5R will vanish in the two limits a — > where s' 12 — > s±2, d 12 — > C12, s'^ — > 0, 
and di — > 1, and a — > oo where s' 12 — > 1, c' 12 — > 0, s'^ — > ci%(— S13), and di — > 513(013) for 
normal (inverted) hierarchy. The question is whether 5R will stay negligible in between as 
s' 12 runs from s±2 to 1, c' 12 from C12 to 0, s', from to C13 (normal) or — S13 (inverted), 
and di from 1 to S13 (normal) or C13 (inverted) as shown in Fig. 18. The dependence of 
the non-zero elements of 5R on (3 = — log e (a/|(5m| 1 |) is shown in Fig. 19. The bumps at 
a ~ Srri2i for both hierarchies, and that at a ~ for the normal hierarchy, happen due 

to the 8'i2 factor competing with the (ft' factor as one of them goes through a resonance 
while the other damps to zero. The heights of the bumps depend on the narrowness of the 
resonances. 

For the case shown in Fig. 19, which was generated with the numbers in Table 1 as 
input, all elements of SR are 0(e 3 ) or smaller for the entire range of a, with the maximum 
value of ~ 0.01 ~ 2e 3 occurring in d 12 s'^ near a ~ ^^31 in the normal hierarchy case. Since 
the size of the third rotation angle we neglected in the Jacobi procedure was 0(s 2 sis), 
Eq. (B.31) is valid to the same order provided S13 = 0(e). 

For smaller values of S13, the resonance at a ~ would have been narrower, and the 
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Figure 20. /3-dependence of the non-zero elements of SR for different values of S13 with normal 
hierarchy, (a) s 13 = 0.03 = 0(s 2 ), (b) s 13 = 0.005 = 0(e 3 ). 



peaks in c' 12 s^ and c' 12 (l — c^) higher. This is illustrated in Fig. 20. In the limit S13 — > +0, 
s'^ and 1 — will become step functions at /? ~ 0, and the maximum height of the peak 
will be 



z l2 {a ~ 5m 



31/ 



0.46 e l = 0.014 , 



(C.2) 



as can be discerned from Eq. (B.17). This is the same as the asymptotic value of ac^/^m^ 
discussed earlier. While this value may not seem particularly large, only a factor of 3/2 
larger than the peak in Fig. 19(a), it is parametrically 0(e 2 ). On the other hand, the third 
rotation angle neglected in the Jacobi procedure was 0(e 2 si3). Thus, using Eq. (B.31) 
would lead to dropping terms that are larger than the ones we keep when S13 = 0(e 2 ) or 
smaller. Also, the sudden change in the accuracy of Eq. (B.31) across a ~ Sm^, as can be 
seen in Fig. 20, will lead to kinks in the resulting oscillation probabilities. 

C.2 Anti-neutrino Case 

The difference between the two sides of Eq. (B.59) is given by 



6R = Rii(6 12 , 0)i?i 3 (0 , 0) - R 13 (cf> , o)R 12 (e 12 , 0) 



'12V 



(1 



-12) 



'i2(l-5y -(l-ci 2 K 
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-s'12^ 



S 'l2 s '<j> 
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It is clear that 5R will vanish in the two limits a — > where s'i 



and c'< 1, and a — > 00 where s'i 



12 



0, c' 



12 



12 — > s 12, C 12 

1, s', -si 3 (c 13 ), and c!. 



C12, s'^ -> 0, 
► ci 3 (si 3 ) for 

normal (inverted) hierarchy. The question is whether 5R will stay negligible in between as 
s' 12 runs from s± 2 to 0, c' 12 from c\% to 1, from to — S13 (normal) or C13 (inverted), and 
from 1 to C13 (normal) or S13 (inverted) as shown in Fig. 21. 
The dependence of the non-zero elements of 5R on (3 = — log e (a/|5m 31 |) is shown in 
Fig. 22, which was generated with the numbers in Table 1 as input. We can see that all 
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Figure 21. /3-dependence of s' 12 , 1 — c' 12 , s'^ and 1 — c'^ for (a) normal and (b) inverted hierarchies. 
The behaviors of s' 12 and 1 — c' 12 are common to both. 
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Figure 22. /3-dependence of the non-zero elements of 5R for (a) normal and (b) inverted hierarchies. 
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Figure 23. /3-dependence of the non-zero elements of 5R for different values of S13 with inverted 
hierarchy, (a) s 13 = 0.03 = 0(s 2 ), (b) s 13 = 0.005 = 0(e 3 ). 



elements of 5R are 0(e 3 ) or smaller for the entire range of a, with the maximum value of 
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~ 0.01 ~ 2e 3 occuring in s' 12 s^ near a ~ in the inverted hierarchy case. Since the size 
of the third rotation angle we neglected in the Jacobi procedure was 0(e 2 si3), Eq. (B.59) 
is valid to the same order provided S13 = O(s). 

For smaller values of S13, the resonance at a ~ would have been narrower, and the 
peaks in s'^s'^ and s' 12 (l — c^) higher. This is illustrated in Fig. 23. In the limit S13 — > +0, 
s'^ and 1 — will become step functions at (3 ~ 0, and the maximum height of the peak 
will be the same as Eq. (C.2), and the asymptotic value of as'^/l^m^l, as can be discerned 
from Eq. (B.46). This is parametrically 0(e 2 ), while the third rotation angle neglected in 
the Jacobi procedure was 0(e 2 si3). Thus, using Eq. (B.59) would lead to dropping terms 
that are larger than the ones we keep when S13 = 0(e 2 ) or smaller. Also, the sudden 
change in the accuracy of Eq. (B.59) across a ~ Sm^, as can be seen in Fig. 23, will lead 
to kinks in the resulting oscillation probabilities. 
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